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Abstract. This paper treats the solvability of a semilinear reaction-diffusion system, 
which incorporates transport (diffusion) and reaction effects emerging from two separated 
spatial scales: x - macro and y - micro. The system's origin connects to the modeling of 
concrete corrosion in sewer concrete pipes. It consists of three partial differential equations 
which are mass-balances of concentrations, as well as, one ordinary differential equation 
tracking the damage-by-corrosion. The system is semilinear, partially dissipative, and 
coupled via the solid-water interface at the microstructure (pore) level. The structure of 
the model equations is obtained in by upscaling of the physical and chemical processes 
taking place within the microstructure of the concrete. Herein we ensure the positiv- 
ity and bounds on concentrations, and then prove the global-in-time existence and 
uniqueness of a suitable class of positive and bounded solutions that are stable with re- 
spect to the two-scale data and model parameters. The main ingredient to prove existence 
include fixed-point arguments and convergent two-scale Galerkin approximations. 

Keywords: Reaction and diffusion in heterogeneous media, two-scale Galerkin ap- 
proximations, parabolic variational inequality, well-posedness 



1 Introduction 

We consider a two-scale (distributed-microstructurej^ or double-porosity-) system mod- 
eling penetration of corrosion in concrete sewer pipes. This kind of models appears in 
a multitude of real-world applications and are therefore of great importance mainly be- 
cause they are able to connect the information from the microscale to the macroscale 



This terminology is very much due to R. E Showalter; see chapter 9 in |10) . 



(e.g. via the boundary of the cell, micro-macro transmission conditions). They are usu- 
ally obtained in the homogenization limit as the scale of the inhomogeneity goes to zero. 
These models provide a way to represent a continuous distribution of cells within a global 
reference geometry. Roughly speaking, to each point x E fl, we assign a representative 
cell Yx. The flow within each cell is described (independently w.r.t what happens at the 
macroscale) by an initial-boundary-value problem. The solution of the problem posed in 
the cell Yx is coupled via the boundary of Y^ to the macroscale. In order to derive such 
models, different choice of microstructures within the global domain can be considered. 
Note that the microstructure Y^ does not necessarily need to be periodically distributed 
in Q. Uniformly periodic (i.e., Y^ = Y) and locally periodic are some of the options. 
Apparently, this sort of models are good approximation of real situations, but involve the 
computation of a very large number of cells problems. 

The model was obtained via formal homogenization using the different scalings in e 
(a scale parameter referring to the micro-geometry) of the diffusion coefficients in locally 
periodic case; details can be seen in [7]. The relevant mathematical questions at this point 
are threefold: 

(Ql) How does the information 'flow' between macro and micro scales? What are the 
correct micro-macro transmission boundary conditions? Is the model well-posed in 
the sense of Hadamard? 

(Q2) What is the long-time behavior of this two-scale reaction-diffusion system? 

(Q3) To which extent is the model computable? 

In this paper, we focus on the first type of questions - the well-posedness of the two- 
scale model - preparing in the same time the playground to tackle the second question. 
For more information on the modeling, analysis and simulation of two-scale scenarios, we 
refer the reader to HTJ H^l HSl [13 EHl EI]- We postpone the study of the long-time 
behavior (and of the inherently occurring memory effects) to a later stage. Also, we will 
investigate elsewhere the computability of our system. However, it is worth mentioning 
now preliminary results in this direction: ^Lbj proves the rate of convergence for a two- 
scale Galerkin scheme, while the authors of [S] we produce the first numerical simulations 
of two-scale effects to our system. Our working techniques are inspired by [121 HZ] . The 
novelty we bring in this paper is two-fold: 

• We are able to apply the two-scale Galerkin procedure to a partly-dissipative R-D 
system, (compare to the case of fuUy-dissipative system treated in [IS1II7]). 

• We are able to circumvent the use of both extra (macro) x— regularity of the micro- 
solutions and if ^—regularity w.r.t. y of the initial data by making use of the 
parabolic variational inequality framework for part of our system. 

The main reason why we want to keep low the regularity assumption on the initial data 
is that, at a later stage, we would like to couple this reaction-diffusion system with the 
actual degradation of the material, namely with partial differential equations governing 
the mechanics of fracturing concrete. 



The paper is organized as follows: Section [2] contains a brief introduction to the 
chemistry of the problem as well as a concise description of the geometry. The two-scale 
model equations are introduced in Section [3] while Section |4] includes the functional setting 
and assumptions as well as the main results of the paper. The proofs of the results are 
given in Section |6] 



2 Description of the problem 
2.1 Chemistry 

Sulfuric attack to concrete structures is one of the most aggressive chemical attacks. 
This happens usually in sewer pipes, since in sewerage there is a lot of hydrogen sulfide. 
Anaerobic bacteria present in the waste flow produce hydrogen sulfide gaseous H2S. From 
the air space of the pipe, H2S enters the air space of the microstructure where it diffuses 
and dissolves in the pore water. Here aerobic bacteria catalyze into sulfuric acid H2SO4. 
As a next step, H2SO4 reacts with calcium carbonate (part of solid matrix, say concrete) 
and eventually destroys locally the pipe by spalling. The model we study here incorporates 
two particular mechanisms: 

• the exchange of -^25' in water and air phase of the microstructure and vise versa [2], 

• production of gypsum as a result of reaction between H2SO4 and calcium carbonate 
at solid- water interface. 

The transfer of H2S between water and air phases is modeled by deviation from Henry's 
law. The production of gypsum (weakened concrete) is modeled via a non-linear reaction 
rate 1]. Here we restrict our attention to a minimal set of chemical reactions mechanisms 
(as suggested in |4j), namely 



' lOH^ + SO^^ + org. matter — > H2S{aq) + AH2O + oxid. matter 
H2Siaq) + 2O2 2H+ + SOf 

H2S{aq) ^ H2S{g) ^ ^ 

^ 2H2O + H+ + SOf + CaCO^ — > CaSOi ■ 2H2O + HCO^ 

We assume that reactions ([T]) do not interfere with the mechanics of the solid part of the 
pores. This is a rather strong assumption: We indicated already that ([T]) can produce 
local ruptures of the solid matrix [18], hence, generally we expect that the macroscopic 
mechanical properties of the piece of concrete will be affected. For more details on the 
involved cement chemistry and practical aspects of acid corrosion, we refer the reader to 
[3] (for a nice enumeration of the involved physicochemical mechanisms), [18] (standard 
textbook on cement chemistry), as well as to [131 HH] cind references cited therein. For a 
mathematical approach of a theme related to the conservation and restoration of historical 
monuments [where the sulfatation reaction ([T]) plays an important role], we refer to the 
work by R. Natalini and co-workers (cf. e.g. [I]). Based on the reactions mechanisms 
single-scale and two-scale ([!]), some models were derived in [T] El El E] mostly relying 



on periodic and locally-periodic homogenization. The model we investigate here has the 
same structure as one of these cases elucidated via homogenization. 

2.2 Geometry 

We consider Q and Y to be connected and bounded domains in R,^ defined on two different 
spatial scales, is a macroscopic domain with Lipschitz continuous boundary F which is 
composed of two smooth disjoint parts: and F^. F is a standard unit cell associated 
with the microstructure within Q with Lipschitz continuous boundary dY and has three 
disjoint components, i.e., Y := YqUYiIJ Y2, where Yo,Yi,Y2 represent the solid matrix, 
the water layer which clings to the pipe wall and air-filled part surrounded by the water 
in the pipe, respectively. All constituent parts of the pore connect neighboring pores to 
one another. Fi, F2 denote the inner boundaries within the cell Y, see Figure [T| that 
is, dYi = Fi U F2 U {dYi fl dY). Fi represents the solid- water interface where the strong 
reaction takes place to destroy the solid matrix and F2 denotes the water-air interface 
where the mass transfer occurs. Fi,F2 are smooth enough surfaces that do not touch 
each other. 




Figure 1: Left: Cross-section of a sewer pipe. Second from the left: A cubic piece from the 
concrete wall zoomed out. This is the scale we refer to as macroscopic. Second from the right: 
Reference pore configuration. Right: Zoomed one end of the cell. 



3 Two-scale model equations 

The two-scale reaction-diffusion system we have in mind consists of the following set of 
partial differential equations coupled with one ordinary differential equation: 

dtW,-Vy{d^VyW,) = -fi{wi) + f2{w2) in(0,T)xl]xn, (2) 

dtW2-Vyid2VyW2) = fiiwi)- f2iw2) in (0, T) X X Fi, (3) 

dtws - V ■ (dsVws) = -« y (^^3 - W2)d-fy in (0, T) x fi, (4) 

dtWi = r]{wi,W4) on (0, T) X fi X Fl. (5) 



The system is equipped with the initial conditions 



Wj{0,x,y) 
W3{0,x) = 



{x,y), J €{1,2} 



in 17 X Yi, 



in fi, 



^4(0, X, y) = wl{x, y) on 11 X Ti, 



(6) 



while the boundary conditions are 



diWyWi ■ u(y) = —r](wi,W4) on (0,T) x x Fi, 
diVyWi ■ v{y) = on (0, T) X X T2 and (0, T) x fi x {dYi n dY), 
d2VyW2 ■ v{y) = on (0, T) X fi X Ti and (0, T) x fi x {dYi n dY), 
d2VyW2 ■ v{y) = a{Hw3 - W2) on (0, T) x fi x 
dsVwa ■ u{x) = on (0, T) x T^, 
W3 = w^ on{0,T)xrD, 



(7) 



where Wi denotes the concentration of H2S0i in (0,T) x fi x Yi, W2 the concentration of 
H2S aqueous species in (0,T) x fi x Fi, W3 the concentration of H2S gaseous species in 
(0,T) X fi and W4 of gypsum concentration on (0,T) x fi x Fi. V without subscript de- 
notes the differentiation w.r.t. macroscopic variable x— , while Vy,divy are the respective 
differential operators w.r.t. the micro-variable y. a denotes the rate of the reaction taking 
place on the interface r2 and H is the Henry's constant. The microscale and macroscale 
are connected together via the right-hand side of and via the micro-macro boundary 
condition ([t])^. The information referring to the air phase Y2 is hidden in w^. The partial 
differential equation for w^, defined on macroscopic scale, is derived by averaging over Y2; 
details are given in [7|. 



4 Functional setting. Assumptions. Main results 

In this section, we enumerate the assumptions on the parameters and initial data needed 
to deal with the analysis of our problem. Furthermore, the definition of the solution to 
the system (|2|-([7]) is discussed and the main results of the paper are given at the end of 
this section. To keep notation simple, we put 

X ■={ze H\n)\z = on Yd}. 



Assumption 4.1 (Al) d, G L~(fi x Yi),i e {1,2} and d^ G L°°(fi) such that 

{di{x,y)^,^) > d^\C,\'^ for d^ > for every ^ e R^, a.e. {x,y) eQxYi audi G {1,2}, 
and {ds{x)C,,^) > (i^l^P for d^ > for every ^ G and a.e. x G fi. 

(A2) ri{a,P) := R{a)Q{P) , where R andQ are locally Lipschitz continuous functions such 
that R' > and Q' < a.e. on R and 

R(a) ■= I ^/ " > Q(g\ f positive, if (3 < (3max, 

\ 0, otherwise, I 0, otherwise, 

where Pmax is a positive constant. Also, we denote by R the primitive of R with 
R{0) = 0, that zs, R{r) = Jq R{Od^ for r G R. 



(A3 ) The functions fi,i G {1,2}, are increasing and locally Lipschitz continuous functions 
with fi{a) = for a < and fi{a) > for a > 0, i E {1,2}. Furthermore, 
TZifi) = 'R-{f2), where TZ^f) denotes the range of the function f . Obviously, for 
Ml, M2 > there exist positive constants M[, M'2 > such that 

fi{M[) = /2(M2), M[ > Ml and M'^ > M2. 

(A4) wio G L2(fi; H\Yi)) n L^{n x Y^), W20 G L^Q; H\Yi)) f] L^{n x Y^), 

W30 G H^{n)nL^{n), W3o-w^{0,-) eX,w^ e L^{0,T;H^{n))nH\0,T;L\n)) 
n L^((0,r) X n) with Vwf ■ z/ = O on (0,r) x Tn, W40 G X Ti). 

Note that in (A4) we define L^(n') := L°^{n') n {u\u > on Q'} for a domain n'. 
Next, we give the definition of a suitable concept of solution to our problem: 

Definition 4.2 We call the multiplet {wi,W2,W3,W4) a solution to the problem if 
(SI) ~ (S5) hold: 

(51) WuW2 G H\0,T;L'^{n x Y^)) n L'^{0,T; L'^{n; H\Yi))) f] L^{{0,T) xQx Y^), 
W3 G H\0,T;L^{n)) D L°^((0,T) x n), W3 - e L~(0,T;X), 

W4 G H\0,T;L\n x Ti)) n L°°((0, T) xQx Ti). 

(52) It holds that 

dtWi{wi — vi)dxdy + / diVyWi ■ \/y{wi — vi)dxdy 

QxYi JqxYi 

+ / Q{w4){R{wi) - R{vi))dxdjy 

< / {-fi{wi) + f2{w2)){wi - vi)dxdy 
JnxYi 

forvi G L'^{n;H\Yi)) with R{vi) G L^{VL x Ti) a.e. on [0,T]. 

(53) It holds that 

I dtW2V2dxdy + / d2VyW2 ■ VyV2dxdy — a I {Hws — W2)v2dxd'yy 

JnxYi JqxYi JnxT2 

= I {fi{wi)-f2{w2))v2dxdy forv2 G L\n-H\Yi)) a.e. on [0,r]. 
JnxYi 

(54) It holds that 



nxYi 
a 



dtWsVsdx + / dsVyWs ■ Vv^dxdy 
Jn 

' / {Hws — W2)vsdxd'yy for V3 G X a.e. on [0,T]. 
Jnxr2 

(S5) (0 holds a.e. on (0,T) xVtxTi. 



Theorem 4.3 (Uniqueness) Assume (Al)-(A4), then there exists at most one solution 
in the sense of Definition \4-S\ 



Proof. For the proof, see Section |6} 

Remark 4.1 Having in view the proof of Theorem \4.3\ and the working techniques in 
Theorem 3, pp. 520-521 in 16] as well as Theorem 4-1 in [141, we expect that the solution 
in the sense of the Definition \4.2\ is stable to the changes with respect to the initial data, 
boundary data, and model parameters. 

Theorem 4.4 (Global existence of solutions to J2])-J7|j Assume (Al)-(A4), then there 
exists a solution {wi,W2,wz,W4) of the problem ([2|-([7^. Moreover, it holds that 

(i) Wi{t),W2{t) > a.e. in fl x Yi, W2,{t) > a.e. in VL and w^{t) > a.e. on VL xVi 
for a.e. te [0,T]. 

(a) Wi{t) < Ml, W2{t) < M2 a.e. in x Yi , Ws{t) < M3 a.e. in and W4{t) < M4 
a.e. onQ xTi for a.e. t G [0,T], where Mi, M2, M3 and M4 are positive constants 
satisfying Mi > ||wiio|U°°(nxyi)' ^2 > \\w2o\\L^{nxYi), 

M3 > max{||M;3o||Loo(f7), llLoo(Qxyi)}, /i(A^i) = /2(M2) and M2 = HM3 and 
M4 = max{(3max, |k40||L°°{Qxri)}- 

In order to prove the existence of a solution, we first solve the following problem 
Pi{g, h) in Lemma 5.1 for given functions g and h: 

dtwi - Vy ■ (diVyWi) =g in {0,T) xQxYi, 

diVyWi ■ v{y) = -hR{wi) on (0,T) xVLxTi, 

diVyWi ■ v{y) = on (0, T) X 1] X T2 and (0, T) x VL x {dYi n dY), 

WiifS) = Wio on fl X Yi. 

Next, for a given function g on (0,T) x Q x Yi, we consider the following problem 
P2(fl') (see Lemma 5.2): 

dtwi - Vy ■ {diVyWi) = g in (0,T) x ^] x Fi, 

di'VyWi ■ z/(y) = —ri{wi,Wi) on (0,T) x Q x Ti, 

diVyWi ■ v{y) = on (0, T) X X T2 and (0, T) x Vt x [dYi n dY), 

dtWi = ri{wi,W4) a.e. on (0, T) x Q x Ti, 

Wi{0) = Wio on Q X Yi and ^4(0) = W40 on Q x Ti. 



As a third step of the proof, we show the existence of a solution of the following 
problem P3(5') for a given function g on (0,T) x i7 x "Ki (see Lemma 5.3): 

dtW2 - Vy ■ {d2VyW2) =g m{0,T) xQxYi, 

dtW3 - V ■ (rfsVws) = y (^^3 - W2)d'jy in (0, T) x Q, 

d2VyW2 ■ z/ = on{0,T) xnxTi and (0, T) x n x {dYi n dV), 

d2VyW2 ■ y = a{Hw3 — W2) on (0,T) x Q x T2, 

dsVws ■ u{x) = on (0, T) x T^, 

W3 = on (0, T) X Tfl, 

^2(0) = W20 on Q X Yi and ^3(0) = W30 on f2. 



5 Auxiliary lemmas 

Lemma 5.1 Assume (Al), (A2), (A4), h G H^{0,T; L'^{n x Li)) n L^((0, T) xQxTi) 
and g G L^((0,T) x x yi). // R is Lipschitz continuous and bounded on R, then 
there exists a solution wi of Pi{g,h) in the following sense: Wi G H^{0,T] L^{Q x Yi)) fl 
L^{0,T;L\n;H\Yi))) satisfying 

/ dtWi{wi — Vijdxdy + / diVyWi ■ Vy{wi — Vi)dxdy 
JnxYi JnxYi 

+ / h{R{wi) — R{vi))dxd'^y 

< / g{wi-vi)dxdy for vi e L'^iil; H\Yi)) a.e. on [0,T], (8) 
JnxYi 

wi(0) = 1^10 on Q X Yi. 

Proof. First, let {Q} be a Schauder basis of L^{Q] H^(Yi)). More precisely, {Q} is an 
orthonormal system of a Hilbert space L'^{Q x Yi) and is a fundamental of L^(f2; i7^(yi)), 
that is, for any z G L^{Q x Yi) we can take a sequence {z^} such that z^ = Ylf=i '^jO 
Zk z m. L"^ [VL] {Yi)) as A; — )■ 00, where G R. Then there exists a sequence {tf^Q} 
such that w^Q := Xlj^i o^ioCj a^id iu"q — )■ i^io in LF' [VL] H"^ {Yi)) as n — )■ cxd. 

Here, we are interested in the finite-dimensional approximations of the function wi 
that are of the form 

X, y) := J] a,"(t)C(x, y) for (t, x, y) G (0, T) x x Fi, (9) 

where the coefficients a", j = 1, 2, . . . , iV„, are determined by the following relations: For 
each n 



nxYi 



(dtw'l{t)^i + diV yw'l{t)V y(i)i)dxdy + 



hR{w'l{t))(t)idxd'^y 



= / g{t)4)idxdy for 0i G span{0 : ^ = 1,2, .., A^„} and t G (0,T], (10) 

JqxYi 

a]iO) = a]^ forj = l,2,...iV„. (11) 

Consider (pi = Q, j = l,2,...Nn, as a test functions in (10). This yields a system of 
ordinary differential equations 

dta]{t) + ^(A,),«r W + F^it, «"(t)) = J,(t) for t G (0, T] and j = 1, 2, . . . , iV„, (12) 

where a"(t) := « (t), • . . , a^„(t)), (Ai)i := /^^^y^ rfiVj,Ci ■ VyQdxdy, 
F^{t,a^) := 4^r,^W^(EiI"iO)Ci^^a;rf7j/ and J,{t) = j^^y^g{t)Qdxdy for t G (0,T]. 
Note that F" is globally Lipschitz continuous due to the assumption of this lemma. 
According to the standard existence theory for ordinary differential equations, there exists 
a unique solution a", j = l,2,..,Nn, satisfying (12) for < t < T and (11). Thus the 
solution w'^ defined in ^ solves (lO)-(ll). 

Next, we show some uniform estimates for approximate solutions with respect to 
n. We take 0i = Wi in (10) to obtain 

dtw'l{t)w'l{t)dxdy+ [ di\Vyw'l{t)\'^dxdy+ [ h{t)R{w1{t))w'l{t)dxdjy 



flxYi 



nxYi 



nxTi 



g{t)w'l{t)dxdy fortG(0,T]. 



Q,xYi 



Since R{r)r > for any r G R, we see that 

1 d 
2dt 



nxYi 



\VyWi{t)\^dxdy 



nxYi 



w'l{t)\^dxdy + d\ 

<\l W)?dxdy+ \ f \w'1{t)\^dxdy fortG(0,r]. 
Applying Gronwall's inequality, we have 

w1{t)\^dxdy + d\ f [ \Vyw'l{t)\'^dxdy <C for t G (0,T], 

Jo JqxYi 



nxYi 



where C is a positive constant independent of n. 

To obtain bounds on the time-derivative, we take 0i = dtWi as test function in (10). 
It is easy to see that 

1 d 



/ 



\dtw'l{t)\'dxdy + 
g{t)dtWi{t)dxdy 



2dt 



di\Vywl{t)ydxdy + 



nxY-i 



h{t)dtR{w"l{t))dxd-iy 



nxYi 



< 



\f \g{t)\^dxdy + \ f \dtw'1{t)\^dxdy for a.e. t G (0, T]. 
2 JnxYi 2 jQy^Yi 



Accordingly, we have 

- I \dtWi\^ dxdy + -— I di\V yWi\^ dxdy + — / hR{wi)dxd'jy 
2 jQxYi 2 dt jQxyi dt j^xri 

^ / |(7|^(ix(i?/ + / dthR{w^)dxd'jy a.e. on G (0,T]. 
By integrating the latter equation, we have 

Iff \dtw'l\^dxdydt + ^ [ \Vyw'l{ti)\^dxdy+ [ h{t)R{w'l{ti))dxdjy 

2 Jo JqxYi 2 jQy^Yi JuxVi 

< \ I di\Vyw'^{Q)\^dxdy+ [ h{0)R{w'l{0))dxdjy 

+ ]- j j \g\^dxdydt+ j j dthR{w'l;)dxd'yydt fortiG(0,T]. 
2 Jo JqxYi Jo JnxFi 



'0 JnxYi Jo JQxF 

Note that 



[ [ \R{w'^)\^dxd-iydt < C [ [ \w1\^dxd-iydt 
Jo JnxFi Jo JnxTi 



T 



< C I I (|Vj,<|2 + \w'l\^)dxdydt. 
JnxYi 



Here, we have used the trace inequality. Hence, we observe that {w"} is bounded in 
H^{0,T; L'^{n X Yi)) and L°°(0,T; L'^{n; H^iYi)). From these estimates we can choose a 
subsequence {nj} of {n} such that w^' — )■ Wi weakly in H^{0,T; x Yi)), weakly* in 
L°°(0,T;L2(fi X Fi)), and weakly* in L°^{0,T; L^{n; H^{Yi))). Also, the above conver- 
gences implies that %vl"{T) — t- wi{T) weakly in LF'{Vt x Yi). 

Now, in order to show that ^ holds let v e ^^(0, T; L'^{VL; H^{Yi))). Obviously, we can 
take a sequence {ffc} such that ffc(t) := Xlj^i "^j (^)O aiid — > f in L^(0, T; if-'^(Fi))) 
as — )■ oo, where d^ e C([0, T]) for z = 1, 2, . . . , and A; = 1, 2, . . .. For each k and i 
with iV„. > mfc from (jTol) it follows that 



T 



- v\)dxdydt + / diVyW^' ■ Vy{w'l' - v\)dxdydt 
JnxYi 



JnxY^ 

^ r 

/ hR{w'l^){w'l^ - v'l)dxd^ydt 

JqxFi 

< [ [ g{wi' - v\)dxdydt. 

Jo JqxYi 



By the lower semi-continuity of the norm and the convex function R we have 

liminf ( [ [ {dtw1^{w1^ - v'l) + diVywl^ ■ Vy{w1^ - vl)}dxdydt 



+ / hR{w1'){w1' - v^)dxd'jydt 

Jo JuxVi 



> 



/ / dtWi{wi — Vi)dxdydt + / diVyWi ■ \/y{wi — Vi)dxdydt 
Jo JnxYi JnxYi 

+ / h{R{wi) — R{vi))dxd'-fydt for each k. 
Jo JnxTi 



Then we show that ([8j) holds for each Vk- Moreover, by letting — )■ oo we obtain the 
conclusion of this lemma. □ 

Next, we solve the problem P2(fi')- 

Lemma 5.2 Assume (Al), (A2), (A4), and g e L'^{{0,T) x Q xYi). If R and Q are 
Lipschitz continuous and bounded on R, then there exists a solution {wi,W4) of P2{g) 
m the following sense: Wi G H^{Q,T; L'^iVt x Yi)) n L°°(0, T; ^^(fi; ifi(Fi))) and G 
H\0,T;L'^{n X Ti)) satisfying 

/ dtWi{wi — Vijdxdy + / diVyWi ■ Vy{wi — Vi)dxdy 
JnxYi JnxYi 



+ / Q{w4){R{wi) - R{vi))dxd^y 
JnxFi 

< I g{wi-vi)dxdy for vi e L^{Q; H\Yi)) a.e. on [0,T], (13) 
JnxYi 

dtWi = ri{wi, Wi) on (0, T) x ^] x T, (14) 
Wi(0) = WiQ on Q X Yi and 104(0) = W40 on Q x Ti. 

Proof Let e V := {z e H\0,T;L'^{Q x Fi)) : z{0) = w^o}- Then, since 
Q{w4) e H\0,T;L'^{n X Fi)) n L^((0,T) x x Fi), Lemma |5TT| implies that the prob- 
lem Pi{g,Q{tV4)) has a solution wi in the sense mentioned in Lemma 5.1, Also, we put 
W4,{t) := Jq ?7(wi(r), W4(r))(ir -|- 1^40 on fi x Fi for t E [0,T]. Accordingly, we can define 
an operator A^ : V ^ V hj At(w4) = w^. 

Now, we show that At is a contraction mapping for sufficiently small T > 0. Let 
wi e H\0,T;L\n X Fi)) and w\ be a solution of Pi(^,Q(w^)) and wl = Ariwl) for 
i = 1,2, and Wi = wl — wf, w^ = wl — w\ and W4 = w\ — w\. 

First, from (Is]) with vi = wl we see that 



/ diWywl ■ Vy(w| — wl)dxdy 
JnxYi 

+ 1 Q{wl){R{w\) - R{wl))dxd-iy 



dtWi{wi — Wi)dxdy + 

nxYi JQxYi 



nxFi 



< 



/ g{wl — wl)dxdy a.e. on [0,T]. 
JnxYi 



Similarly, we have 

/ dtw\{w\ — w\)dxdy + / diVyVol ■ Vj,(wi — w\)dxdy 
JnxYi JnxYi 

+ [ Q{wl)iR{wl) - R{wl))dxd-fy 

< / g{wl — wl)dxdy a.e. on[0,T]. 
JnxYi 

By adding these inequalities, for any £ > we obtain 

\wi\^dxdy + / di\\/yWi\'^dxdy 

JqxYi 

< - [ {Q{wl) - Q{wl)){R{wl) - R{wl))dxd-fy (15) 

< Ce {wil'^dxd'jy + e / {wil'^dxd'jy 

Jnxri Jnxri 

< Ce Iw^l'^dxd'jy + CyiE / (iWyWil"^ + \wi\'^)dxdy a.e. on[0,r], 

Jnxri JqxYi 

where Cy^ is a positive constant depending only on Yi. Here, by taking e > with 
Cy^e = ^d1 and using Gronwall's inequality we see that 

If d^ 

' \wi{t)\'^dxdy ' ^ 



1 d 
2dt 



nlVyWil'^dxdydr 
* /" [ \w4'^dxd-fydT forte[0,T] 

Jo JuxFi 

Next, on account of the definition of it is easy to see that 

\w4{t)\'^dxd'yy 



2 .'nxYi 2 
< e^'* / / \w4'^dxd-fydT forte[0,T]. (16) 

QxFi 



1 d 
2dt 



< \ I (|r7(ti;|(t),^l(t))-r7(u72(t),^2^t))|2 + |^^(t)|2)^3,^^^ 



2 



< cl {\wi{t)\'^ + \wi{t)\^ + \wi{t)\^)dxd-iy for a.e. t e [0,T]. 

JnxVi 



Gronwall's inequahty, viewed in the context of (16), implies that 

\wis^{t)\'^dxd'~^y 



t 

Ct 



n\wi\^dxd'^ydT + I I {\V yWi\^ + \wi\^)dxdydT) 



< 



Ce^* /" [ Iml'^ dxd-^ydr fortG[0,r]. 
Jo JnxFi 



Hence, we obtain 



\dt'Wi\\L^o,T;L^{nxri)) < Wvi'^l,^!) -v{wl,'^l)\\L^o,T;L^{nxri)) 

< C{\\Wi\\L2(^0,T;L2(n;m{Yi)) + || ""^4 1| L2(0,T;L2{Qxri)) ; 

< C||w4||i2(o,T;L2(Cxri)) 

< CT^/^||9(U'4||i2(o,T;L2(nxri))5 



and 



||At(u^1) - Ar(ti'4)||Hi(0,T;L2(Qxri)) < ||'y^4||L2(o,T;L2(Cxri)) + || <9tW^4 || L2 (0,T;L2(f7xri)) 

< Cr^/^||w4||//i(o,T;L2(Qxri))- 

This concludes that there exists < Tq < T such that A^o is a contraction mapping. 
Here, we note that the choice of Tq is independent of initial values. Therefore, by applying 
Banach's fixed point theorem we have proved this lemma. □ 



As the third step of the proof of Theorem 4.4, we solve P3{g)- 



Lemma 5.3 Assume (Al), (A2), (A4), and g e L^((0,T) x Q xYi). Then there exists 
a pair {w2,W3) such that W2 G H^{0,T;L'^{n x Yi)) n L^{0,T; L^{n; H\Yi))), G 
H\0, T- L'^{n)) n L°°(0, T; H\n)), 

/ dtW2V2dxdy + / d2VyW2 ■ VyV2dxdy — a I {Hws — W2)v2dxd'-)y 
JqxYi JqxYi Jnxr2 

= / gv2dxdy for V2 e L'^{Q; H\Yi)) a.e. on [0,T]. (17) 



and (S4)- 



Proof. Let {Q} be the same set as in the proof of Lemma 5.1 and {yUj} be an orthonormal 
system of the Hilbert space L'^{Q) and a fundamental of X. Then we can take sequences 
{w20n} and {Wson} such that w'!^q := J2f=iPjoCj^ ^30 '■= ^f=iPjof^j^ "^20 ^ ^20 in 
L^{Q; H\Yi)) and W^^ W30 - w^?(0) in X as n 00. 

We approximate W2 and W3 := W3 — by functions W2 and of the forms 

w-{t) = J2Pjit)(v W,^it) = J2ljm forn, (18) 

where the coefficients Pj' and 7", j = 1, 2, . . . , A^„ are determined the following relations: 
For each n , we have: 

dtW24>2dxdy + / d2'VyW2 ■ Vy(j)2dxdy — a I {HW^ — W2)4>2dxd'yy 
g(j)2dxdy + a Hw^ (f)2dxd'yy (19) 

nxYi JnxT2 

for 02 G span{Ci : i = l,..,Nn},t G (0,T], 

= f3%, 



Yi Jn JQxr2 

[ {dtw^ - Vdz{Vw^))(t)zdxdy -a [ Hw^ 



(p^dxd'^y 



(20) 



for 03 G span{/ii : i = 1,2, A^„},t G (0,T], 
77(0) =7j^ for j = l,2,..,iV„. 



Consider 02 = Cj and 03 = /i^, j = 1,2, ..,A^„, as a test functions in (19) and (20) 
respectively, these yield a system of ordinary differential equations 



N„. 



9tfi?it) + Y.^B,),/3:{t) + ^(5,),7r(t) = Mt) for t e (0, T] and j = 1, 2 

i=l i=l 

Nn N„ 

9a7it) + J2^C,),j:{t) + 5^(a),/3r(t) = Mt) for t e (0,T] and j = 1,2 



, . . . , iVn, 



i=l 



i=l 



where (5,)^ := J^^y^ rfsV^O ■ VyQdxdy, {Bi)j := J^^^^ fiiQdxd-fy, 

(COi '■= In d-i^^-i ■ V^ijdx + aH J^^^^ fiifijdxd-fy, {Ci)j := -a J^^^^ fijQdxd-fy, 

J^jit) := f^^y^g{t)Cjdxdy, .hj{t) := f^^y^{dtw^ -V{Vw^))Qdxdy+aH f^^^^w^Qdxd-fy 

for t e (o,r]. 

Clearly, this linear system of ordinary differential equations has a solution /3" and 7^". 



Thus the solutions W2 and defined in (18) solve (19) and (20), respectively. 

Next, we shall obtain some uniform estimates for W2 and W3 • We take 02 
03 = in (19) and (20), respectively, to have 



W2 and 



J4 / \w^{t)\'^dxdy+ f d2\VyW^{t)\'^dxdy + j I |>^2vvi— 



nxYi 



\w2it)\'^dxd'yy 



< ^H\T2\ [ \W^{t)\'dx + l [ mi'dxdy 
Jn ^ JnxYi 

+ - I \w'^{t)\'^dxdy + aH^ I \wf{t)\^dxdy for a.e. t G [0, T], 
2 JnxYi JnxV2 



1 d 
2dt 



\W!^{t)\^dx+ / d3\VW^{t)\^dx + aH\T2\ / \W!^{t)\^dxd^, 



< 



nxr2 



\w^{t)\^dxd^y + {a^\T2\ + -) / \W^{t)\^dx 

2 Jn 



- [ \g^{t)\^dx for a.e. t G [0,r], 
2 Jn 



where |r2| := J-p^ d'jy, and := dtw^ — V(Vwf') + aH\T2\w^ . By adding these inequal- 
ities, we get 

~ I \w-2{t)?dxdy + dl [ \VyW^it)\'dxdy 

I \wm'dx+di I \VWm'dx 
< («'|r2| + ^ + f ir^l) / \wm'dx + \ I umxdy 

+ 7 / \w:^it)\^dxdy+l [ \g^{t)\^dx for a.e. t G [0, T]. 
4 JnxYi 2 Jq 

Consequently, Gronwall's inequality implies that for some positive constant C 



I \w^{t)\'^dxdy + I \ W^{t)\^dx <C for t G [0,T] and n, (21) 
JqxYi Jn 

[ [ \VyW^\^dxdydt+ [ [ \VW^{t)\'^dxdt < C for n. (22) 
Jo JnxYi Jo Jn 



'nxYi 

To obtain a uniform estimate for the time derivative by taking 02 = dtW2 in (19) we 



observe that 



/ \dtW^{t)\'^dxdy + ~ [ d2\VyW^{t)\^dxdy 
JnxYi ^ JnxYi 



{H{W^{t) + (t)) - w'i{t))dtw^{t)dxd% 



y 



nxr2 



/ g{t)dtw'^{t)dxdy for a.e. t G [0, T]. (23) 
JnxYi 



Here, we denote the third term in the left hand side of (23) by J{t) and see that 



At) ■■= -^J H{w;\t) + w^{t))w-{t)dxd^y + ~ j Wim^dxd^y 

+H [ dt{W^it)+wl^{t))w^{t)dxd-fy for a.e. t G [0,T]. 

Then we have 

\ f \ d f \ d f 

- \dtw^{t)\^dxdy + -- d2\VyW^{t)\^dxdy+-- \w^{t)\^dxd^y 

<\l \g{t)\'dxdy + ^ [ H{W^{t)+w^{t))w^{t)dxd^y 

+H\T2\'^^ [ \dtiW^{t) + w!^m\\w^{t)\\L^r2)dx for a.e. t G [0,T]. 
Jn 



Similarly, by taking 03 = dtW^ in (20) we have 

m,-{t)\'dx+~Jjs\'^w^it)\'dx 



-a / {H{W^{t) + wf{t)) - w^{t))dtW^{t)dxd-fy - / g^{t)dtW^{t)dx 
Jnxr2 Jq 



< - 



aH\T2\f- [ \W^{t)\^dx + 2a^H^\T2\^ [ \w^{t)\^dx 

2a^\T2\ [ \wJ^it)\'^dxd-fy + l [ \dtW^{t)\^dx+ [ \g^{t)\^dx ior a.e. t e [0,T]. 



+ 

From these inequalities, it follows that 



\dtw:i{t)\'dxdy + -- d2\Vywl^{t)Ydxdy+ -— \w'i{t)\'dxd-iy 

+ \ I \dtW-{t)\'dx+\j^ I d,\VWm'dx + aH\T2\j^ j \Wm'dx 

<\f \g{t)\'dxdy + ^ f H{W^{t)+w!^{t))w^{t)dxd^y 
^ JnxYi Jnxr2 

+ {2\T2\+2a^\T2\) [ \w:^{t)\^dxd-fy+[\g^{t)\^dx for a.e. t G [0, T]. 



'Cxr2 

Here, we use Gronwall's inequality, again, and have 

-ti 



If [ \dtw^\^dxdydt+l [ d2\VyW^{h)\^dxdy+^ I \w'^{h)\''dxdjy 
^ Jo JnxYi ^ JnxYi 



nxr2 



< e 



I r I \dtW^\'dxdt + \ [ d,\VW^{ti)\'dx + aH\r2\ [ m{ti)\'dx 
° Jo Jn ^ Jn Jnxr2 

r ( [ \g\'dxdy+ [ Ig^'l'dx^dt 
Jo \JnxYi Jn / 

+e^*^ ['e-^\^ [ H{W^ + wf)w^dxdjy)dt 
Jo dt Jnxr2 

\g\^dxdy+ [ \g^\'^dx)dt + I H(W^{ti) + wf{ti))w:^{ti)dxd'yy 
\ Jn JnxT2 



10 JnxYj 

+e^*^| [ H{W^{0)+wl^{0))w^{0)dxd-fy\ 
JnxT2 

+e^*^ f ' [ H{W^ + wi^)w^dxd^ydt for h G [0, T]. 
Jo JnxV2 

This inequality together with (21) and (22) leads to 



\VyW':i{t)\^dxdy + f \\/WJ}{t)\^dx <C for t e [0,T] and n, (24) 
Yi Jn 

[ [ \dtw^\^dxdydt+ I [ \dtW^\^dxdt < C hi n. (25) 
Jo JnxYi Jo Jn 



By (21) ~ ( |25) ) there exists a subsequence {ui} such that W2 weakly in 

H^iO,T; X Yi)), weakly* in L°°(0,T; L'^{n; H\Yi)) and W^' -> IVg weakly in 
H\0,T;L\n)), weakly* in L°°(0, T; //^(fi)) as 
L^((0,T) X i7 X as i — i- 00. Here, we put - 



i — > 00. Clearly, 



Wr, 



— )■ W2 weakly in 



Since the problem P3(5') is linear, similarly to the last part of the proof of Lemma 5.1 



we can show (17) and (S4). 



□ 



6 Proof of our main results 

First, we consider our problem (|2])-([7]) in the case when /i, /2, R and Q are Lipschitz 
continuous and bounded on R. 

Proposition 6.1 If (A1)-(A4) hold and fi, f2, R and Q are Lipschitz continuous and 
bounded on R, then there exists one and only one multiplet (^1,^2,^3,^4) satisfying 

Wi,W2 e H\0,T;L^{n x Yi)) f] 1^=^(0, T; L^Q; H\Yi))), 
(S') { W3 e H\0,T-L\n)),W3-w^ e L°°{{0,T-X),W4 e H\0,T-L\n xTi)), 
(S2) holds for any Vi e L^Q; H\Yi)), and (S3), (S4) and (S5) hold. 



Proof. Let (^1,^2) ^ -^^((0,7") x Q x Yi)^. Then, by Lemmas 5.2 and 5.3, there 
exist solutions (^1,^4) of P2(-/i(wi) + /2(w2)) and_(w;2,W3) of P3(/i(w;i) - /2(w2)), 
respectively. Accordingly, we can define an operator At from L^((0,T) x f2 x Yi)"^ into 
itself. From now on, we show that Ay is contraction for small T. To do so, let {w\,W2) ^ 
L2((0,T) X n X {w\,wi) and (wlwl) be solutions of F2{-fi{w{) + f2{wi)) and 

P3(/i(^i) — f2{w\)), respectively, for i = 1,2, and put wi = w\ — wl, W2 = wl — wl, 
Wj = w] - w], j = 1,2,3,4. 



Similarly to (|15), we see that 

/ \wi\'^dxdy + / dil'VyWil'^dxdy 



1 d 
2dt 



< - / iQ{wl) - Q{wl)){R{wl) - R{w\))dxd^y 
ihi^i) - fi{wl))widxdy + / U^iwl) - f2{wl))widxdy 

< I \V yWi\^dxdy + C f \wi\'^dxdy + C f {w^l'^dxd'fy 

2 JqxYi JnxYi JuxFi 



+C / {\wi\'^ + \w2\'^)dxdy a.e. on[0,T]. 

JqxYi 



Next, we test (17) by W2- Consequently, by elementary calculations, we obtain 
1 d 



2dt 



\w2\^dxdy + d2 / \V yW2\^ dxdy + a I \w2\^dxd'^y 
nxYi JnxYi JQxr2 



< / {fi{wl)-fi{wl))w2dxdy- {fi{wl) - fi{wl))w2dxdy 



+a / Hw3W2dxd'yy 
'nxr2 

< C I {\wi\ + \w2\)\w2\dxdy + 



I \w2\''dxd^y+'^H^\T2\ I Iw^l^'dx 



and 



1 d 
2dt 



/ \w2\'^dxdy + d2 / \'VyW2\'^dxdy + — / \w2\'^dxd'yy 
JnxYi JnxYi 2 jQxr2 

C I {\wi\^ + \W2\^ + \w2\^)dxdy + C ! I^g^cia; a.e. on[0,T]. 



< 

'nxYi 

It follows form (S4) that 
1 d 
2~dt 
a 



< 
- 4 



\w'^\^dx + d^ I \Vw'i\^dx + aH I \w'j,\^dxd'^y 
n Jn Jnxr2 

/ \w2\'^dxd'yy + a\T2\ / Itys^dx a.e. on[0,T] 
Jnxr2 Jn 



Moreover, by using the trace inequality and (14), we see that for £ > we can write 
1 d 



2dt 



1^4 1 dx 
nxFi 



< / \r]{wl,wl) - r]{wl,wl)\\w4\dxd-fy 
JnxFi 

< C / (|wi||w4| + \wi\'^)dxdjy 

JnxFi 



< Cy^s / {iVyWil"^ + \wi\'^)dxdy + C / {w^l'^dxd'yy a.e. on[0,T]. 
JnxYi JnxFi 

dP 

Here, we take e with Cy^e = y and add the above inequalities. Then it holds that 

Ty-j; [ \wi\'^dxdy + ^ f \VyWi\'^dxdy 
^ JnxYx 4 JnxYi 

+ / \w2fdxdy^£2 [ \VyW2\'^dxdy + J f \w2\'^dxd'yy 
2"iJqxYi JnxYi JnxT2 

+ 7.-T: I \w2,\^dx + d\ f iWwsl'^dx + aH f \w3fdxd'yy + ]-^ f {w^l'^dx 
^dt Jq^ Jq JnxT2 -^dt Jo^y^Yi 

rr. |2 , I - |2\ T_ 7 , , ^ / /I ., |2 , I .. |2\ 



< C {\wi\' + \w2ndxdy + C {\wiY + \w2V)dxdy 
JnxYi JnxYi 

+C I \w'i\^dx + C I \wA\^dxdy a.e. on [0,T]. 
Jn JnxTi 



Hence, Gronwall's inequality implies that 



{\wi{t)\^ + \w2it)f)dxdy + / \w3{t)\^dx+ / \w4{t)\^dxd'yy 
QxYi Jn JnxFi 



/ / {\wi\'^ + \w2\'^)dxdydT for t e [0,T]. 
Jo JnxYi 



This leads to 

\\ATiwl,wl) - ATiwj,wl)\\Lmo,T)xnxYi) 

< ll^l||L2((o,T)xCxyi) + \\w2\\L'2{{0,T)xnxYi) 

< Ce^'T'/mwl,wl) - iwlwl)\\L2ao,T)xnxn) 

Therefore, there exists a positive number Tq such that Aj^q is a contraction mapping for 
< Tq < T. Since the choice of Tq is independent of initial values, by Banach's fixed 
point theorem we conclude that the problem ([2])-([7]) has a solution in the sense of (S'). □ 

Proof of Theorem \4.4\ First, for m > we define /j^, = 1, 2, Rm and Qm by 

„ / \ _ / fiijn) for r > m, t3 ( \ _ j Rijn) for r > m, 

fim{r) := I ^^^^^ otherwise, "^""^ \ /?(r) otherwise. 

Q{m) for r > m, 
Qm{r) ■■= { Q{r) for |r| < m, 
Q(—m) for r < — m. 



Then, for each m > by Proposition 6.1 the problem ([2])-([7]) with fi = fim, /2 = /2 



I2m, 

R = Rm and Q = Qm has a solution {wim, W2m, w^m, w^m) in the sense of (S'). 
Now, for each m we shall prove 

(i) wim,W2m(t) > a.c. ou (0, T) X Q X Yi, W3rn > a.c. on (0,T) X Q and > 
a.e. on (0,T) x x Ti. 

In order to prove (i) we test (S2) by wim + w^T^, where 0~ := — min{0, 0} with 
(t>'^(t>~ = 0. Then we see that 



1 d 
2dt 



\wim\ dxdy + / di\VyWim\ dxdy 
nxYi JnxYi 



+ / Qmiw4){RmiWlm) - Rni{W\m + W^m))'^^<^ly 

JuxFi 

< / Uim{wim) - f2m{w2m))wimdxdy a.e. on[0,T], 



cxy. 



where Rm is the primitive of Rm with -Rm(O) = 0. Note that Rm{wim)—Rm{wim+Wij 
and (/im(wim) - f2rn{w2ra))wim < 0, siucc f2m > OU R. Clearly, 



1 d 
2dt 



/ \wim\'^d,xdy + / (ii|Vj/Wi„|^(ix(i?/ < a.e. on [0,T] 
JnxYi JnxYi 



so that wim > a.e. on (0, T) x Q x Yi. 

Next, because — [w3m]~ £ X, we can test (S3) by — W2m,^ and (S4) by —w^^ to obtain 



2dt 



- 2 



< - / iflmiWlm) - f2miw2m))w2mdxdy -a Hwzm'W2m dxd^y 

JnxYi Jnxr2 

< 77 / \w2m\'^dxd'yy + ^H'^\r2\ / \w^J^dx a.e. on[0,T], (26) 



and 

1 d 
2dt 



a / {Hw3m - W2m)w^^dxd'yy 
'nxr2 



/ \'^3m\'^dx + dl / \\/w^J^dx 
Jn Jn 

< a \w2m\\f^3m\dxd'yy a.e. on [0,T]. (27) 

JnxVo. 



Adding (|26j) and (27) and then applying Young's inequahty, we get 
1 d 



2dt 



/ \w2j'^dxdy+ / \w3j'^dx) + dl / I VyWa^l^rfxrfy + rfg / {Vw^J'^dx 
JnxYi Jn JnxYi Jn 

< (^if^|r2| + a;|r2|) / \w^jj^\'^dx a.e. on[0,r]. 
^ Jn 

The apphcation of Gronwall's inequahty and the positivity of initial data give W2m > 
a.e. on (0,T) x fl x Yi and wsm > a.e. on (0,T) x Q. 
Since 77 > 0, it is easy to see that 

1 d 



2dt 



/ \w^.i^\dxd'yy < a.e. on [0,T]. 
JnxVi 



Hence, we see that w^rn > a.e. on (0,T) x 1] x Fi. Thus (i) is true. 

Next, we shall show upper bounds of solutions as follows: To do so, by (Al) we can 
take Ml and M2 such that 

Ml > \\wio\\L°°(nxYi), M2 ^ ^^^{\\'l^2o\\L°°{nxYi), H\\w3o\\l^(^q), H\\w!^\\L^(^Qy^Yi),} 

and /i(Mi) = /2(M2). Also, we put M3 = M4 = max{f3max,\\^^o\\L^{nxri)} and 
Mo = max{Mi, M2, M3, M4}. Then it holds: 

(ii) For any m > Mq we have Wim{t) < Mi, W2m{t) < ^2 a.e. in x Fi, w^mit) < M3 
a.e. in Q and ti'4m(^) < M4 a.e. on i7 x Fi for a.e. t G [0,T]. 

In fact, let m > Mq and consider wim — {wim — Mi)+, {w2m — M2)'*" and {wsm — M^)^ 
as test functions in (S2) ~ (S4). Then we observe that 



1 d 
2Jt 



/ \{wirn- Mi) + \^dxdy + dl / \Wy{wim- Mi)+\'^dxdy 
JnxYx JnxYx 



(Fl JnxYx 
+ / Qm{Wim){Rm{Wim) - Rm{wvm " (^Im - Mi) + ))(iX(i7y 



f^xr 



< / {-hrn{Wim) + t2m{w2rn)){Wirn- MiYdxdy, {2i 

JnxYi 



1 d 

2Jt 



JnxYi JnxYi 

< I Ulm{Wim) - f2m{w2m)){w2m- M2ydxdy 

'nxYi 

-to / {Hw3m - W2rn) (^2™ - M2)^dxd'yy, 



(29) 



1 d 
2dt 



[ \{Wsm-Msy\^dx + dl [ \V{Wsm-Msy\^)dx 

Jn Jn 
< -a {Hw3„, - W2m){w3m - M3)^dxd-fy a.e. on [0,T]. 



(30) 



Here, we note that Rm{wim) — Rm{wim — {wim — ^i)^) > 0. Adding (28)-(30), we get 



1 d 
2di 



{\{wim- Miy\' +\{w2m- M2rndxdy+ / \{w3m- M^yi'dx 

nxYi Jn 



+ f {dl\V{wi^ - MO+r + rf°|V(ti;2„ - M2)+|2) dxdy + d° / \V{w3m - M^^l^dx 



'nxYi 

< I {-flm{Wirn) + f2m{w2rn)){{Wirn- MiY - {W2rn- M2Y)dxdy 

'nxYi 



(31) 



+a / {{Hw3m - W2m){w2m " ^2)+ + {W2m " Hw3m){w3m - M3y)dxd-fy a.e. On [0, T]. 

JnxT2 



We estimate the first term on the r.h.s of (31) by making use of fim{Mi) = /2m(^2) and 
the Lipschitz continuity of i = 1, 2, as follows: We have 



[-flmiWim) + flrn{Mi) - f2m{M2) + /2m (W2m) ) (^Im " M^^dxdy 



nxYi 



+ I ifmliWlm) - /im(Mi) + f2m{M2) " /2 m (^^^2) ) (w^2m - M2Y dxdy 

'nxYi 

< I {f2m{w2m) - f2m{M2)){Wim- Mi^dxdy 

'nxYi 



+ I iflmiWlm) - flm{Mi)){w2m - M2Y dxdy 

JnxYi 

< c[ i\{w2rn-M2)^\'' + \{w^m-M^yf)dxdy a.e. on[0,T]. 
JnxYi 



We estimate the second term on the r.h.s in (31) as follows: 

a / (Hwsm - HM3 + M2- W2m){w2m - M2)'^ dxd^y 
+a / {W2m - M2 + H{M3 - W3m)){w3m " M^^dxd'Jy 

< aH {wm3- Ms){w2m- M2)^dxd'yy- a \{w2m - M2)^\'^dxd'yy 

Jnxr2 JnxV2 

+a / {w2m - M2){wi,m - M'i)^ dxd^y - aH I \{wzm - M'i)'^\^dxd^y 
JnxV2 JnxV2 

< {aH"^ + a) / \{wsm — Ms)~^\'^dxd'yy a.e. on[0,T]. 

ioxr2 



Now, (31) becomes 



1 d 
2dt 



{\{Wim-Mi) + \^ + \{w2m-M2) + \^)dxdy+ / |(t/73m-M3) + |2rfx 

nxYi Jn 



+ [ {d^,\ViWim - Mi)+|2 + dl\V{w2m ' M2) + |2) dxdy 
+dl f \V{Wsm-Msy\^dx 

Jn 



< C i\iw2m-M2r\'+\iWvm-M,)^\')dxdy 

'nxYi 



+C I \{w3m - Ms^l^dx a.e. on [0,T]. 
Applying Gronwall's inequality, we get 

i\iwimit)-Mi)+\^ + \{w2mit)-M2)'-\^)dxdy+ I \iw3^{t) - M^^ \^ dx < for t > 0. 



nxYi 

u 

Hence, Wim < Mi, W2m < M2 a.e. in Q x Yi and < M3 a.e. in Q for t G (0, T). 

To show that 104^ is bounded on i7 x Fi, we test ^ with (w4m — ^4)''' and using (A2) 
leads to 

2dt / l(^4m - M4) + |^(ix(i7y 
= J Rm{wi)Qm{w4){w4m - M4)^dxd'jy < a.e. on[0,T]. 

This shows that w^m < M4 a.e on (0,T) x Q x Ti. Thus we have (ii). 

Accordingly, by (i) and (ii) {wim,W2m,'W3m,'W4m) satisfies the conditions (SI) ~ (S4) 
for m > Mq. Thus we have proved this theorem. □ 



Proof of Theorem \4.3 . Let 
(SI) ~ (S4). Since all Wij, i = 
a solution of Q-Q with fi = 
constant m. Then Proposition 6.1 guarantees the uniqueness. This proves the conclusion 
of Theorem 14.31 □ 



[wij, W2j, wsj, W4j), j = 1, 2, be solutions ([2])-([7]) satisfying 
1,2,3,4, j = 1,2, are bounded, {wij,W2j,W3j,W4j) is also 
fim, f2 = f2m, R = Rm and Q = Qm for souie positive 
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